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We consider the dynamics of a multi-component scalar field on super-horizon scales in the 
context of inflationary cosmology. We present a method to solve the perturbation equations 
on super-horizon scales, i.e., in the long wavelength limit, by using only the knowledge of 
spatially homogeneous background solutions. In doing so, we clarify the relation between the 
^-j ' perturbation equations in the long wavelength limit and the background equations. Then 

as a natural extension of our formalism, we provide a strategy to study super-horizon scale 
perturbations beyond the standard linear perturbation theory. Namely we reformulate our 
method so as to take into account the nonlinear dynamics of the scalar field. 



61. Introduction 

> 

j~ — ■ Study of multi-scalar inflation is a current topic in cosmology, motivated mainly by the_fact that 

supergravity theories suggest the existence of many flat directions in the scalar field potentially . Obser- 
vationally we will be able to determine the accurate spectrum of primordial density fluctuations in the 
near future, e.g., by the next-generation MAP and PLANCK satellites!!]) . At that time, to test a variety 
of models of inflation, we will need a systematic method to evaluate the spectrum for a wide class of 
inflaton models. 

In this paper, we consider perturbations of a multi-component scalar field in the long wavelength 
q-( limit. That is, we consider modes whose wavelengths exceed the Hubble horizon scale. In the case 

of single-scalar inflation, the evolution of a perturbation in the long wavelength limit is very simple. 
j^Tjl The behavior of the adiabatic growing mode is specified by 1Z C = constant, where 1Z C is the spatial 

curvature perturbation of the comoving hypersurface. In the case of multi-scalar inflation, the same is 
not true. Even if the wavelength of a perturbation exceeds the horizon scale, 1Z C changes in time. The 
evolution of super- horizon scale perturbations foe, multi-scalar inflation has been investigated in some 
particular models that allow analytical treatments^' . However, since the feasibility of analysis of a model 
does not imply the viability of the model, it is necessary to develop a method that has a wide range of 
applicability. In this context, a rather general framework to study this issue has been given by Sasaki and 
Stewarttf under the assumption of the slow rolling evolution of a multi-component scalar field. However, 
it may well happen that some components of the scalar field do not satisfy the slow rolling condition 
during inflation. Furthermore, the slow rolling condition will be eventually violated toward the end of 
inflation when the reheating (or preheating) commences. Hence it is much more desirable to exclude the 
slow rolling assumption. In this connection, Taruya and Nambu have recently discussed a method to 
obtain the general solution for long wavelength perturbationslP . Unfortunately, however, they have not 
clarified several delicate issues associated with super-horizon scale perturbations that are characteristic 
of general relativity. Here we develop a general framework to study super-horizon scale perturbations 
without assuming slow rolling. Note that, although we use the terminology of 'multi-scalar inflation', 
our framework will be valid for any stage of the universe as long as the energy momentum tensor is 
dominated by a multi-component scalar field. Then we extend our formalism so as to include the effect 
of nonlinearity of the scalar field potential. 

Suppose we are given a set of equations that governs the evolution of a system. If we know a complete 
set of solutions to this set of equations, we have no difficulty in constructing a general solution of the 
perturbation around a fixed background solution. By definition, a complete set of solutions contains a 
sufficient number of constants of motion that parametrize the different solutions. If we take a derivative 
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with respect to one of such constants of motion, a solution to the perturbation equations can be obtained. 
The same is true for the cosmological perturbations. If we could obtain a full set of solutions to the 
Einstein equations coupled with the equations of matter fields, it would be trivially easy to spell out the 
perturbation around a fixed background universe. But it is hopeless to expect so. Nevertheless, if we 
restrict our attention to spatially homogeneous and isotropic universes, we may be able to obtain a full 
set of solutions. Then one may expect that the perturbation in the long wavelength limit is obtained 
from the knowledge of this restricted class of solutions. In fact, although restricted by the slow rolling 
assumption, the formula for the curvature perturbation on the comoving hypersurface derived in Ref. ^) 
is a good example of such a case. 

Let us consider an n-component scalar field whose action is given by 

Smatter = ~J • + U ((/>)) , (IT) 

where we have used vector notation to represent the n-component scalar field. That is, 4> ■ 4> — h pq (j) p (j) q , 
where h pq is the metric in the scalar field space. For simplicity we assume h pq = S pq in the discussions 
below, but the generalization to a non-trivial h pq is straightforward, which we shall discuss later. The 
equations for the spatially flat, homogeneous and isotropic configuration are given by 

4> p + 3Hcj) p + C/I p = 0, (1-2) 
+ (!■«> 

where dot means a derivative with respect to the cosmological time t. We adopt the units 8irG = c = 1, 
H := a/ a is the Hubble parameter with a being the cosmic scale factor, and Jj\ p = h pq dU / 'dcj> q . A solution 
to these equations gives a background solution in the context of cosmological perturbation theory. Hence 
we refer to these equations as the background equations. 

Substituting Eq. (1-3) into Eq. (1-2), we have n coupled second order differential equations. Now, 
suppose that we know a complete set of solutions for these equations. The complete set of solutions will 
contain 2n parameters, A Q (a — 1, 2, • • • , 2n) which distinguish the various solutions. One of them, say A 1 , 
just represents the change in the origin of the time coordinate, to. So let us denote the general solution 
by 4>{t + A , A a ), where a runs from 2 to In. Then a naive anticipation is that we would obtain the 2n 
solutions for the perturbation in the long wavelength limit by d<p(t + A , X a )/d\ a , or more concisely by 
dct>(t,\ a )/d\ a where A 1 = t. 

However, things are not so simple. First of all, it is not guaranteed if we can restrict our consideration 
to spatially homogeneous and isotropic configurations from the beginning to obtain the perturbation in 
the long wavelength limit. Furthermore, it is well known that, when dealing with cosmological pertur- 
bations particularly on super-horizon scales, it is essential to make a clear statement about the choice 
of gauge before one talks about the behavior of the perturbation. In the above discussion, it was not 
clarified in what gauge the solutions were given, even if they actually would describe the long wavelength 
perturbation. One of the purposes of this paper is to clarify these points. 

The paper is organized as follows. In §2, we consider the linear perturbation equations in the long 
wavelength limit and clarify its relation to the background equations. In §§2.1, we introduce the basic 
notation and list the basic equations for the perturbation of a multi-component scalar field. In §§2.2, 
we discuss the long wavelength perturbation in the B = Hl = gauge, where B is the shift vector 
perturbation and Hl is the trace part of the spatial metric perturbation (see Eq. (|2T| ) below). In this 
gauge we find the e- folding number N := log a is unperturbed even under the presence of the perturbation, 
hence it can be regarded as the natural time coordinate instead of the cosmological time, t. Then we show 
that there exists a simple relation between the perturbation equations in the long wavelength limit and 
the derivative of the background equations with respect to the parameters {A a } := {N, X a }. This result 
indicates there is a unique correspondence between the choice of the time coordinate for the background 
equations and that of gauge for the perturbation equations. To strengthen our assertion, in Appendix we 
consider the perturbation equations in a couple of other choices of time coordinates and show that this 
is indeed the case. Then using this relation, we give a formula to calculate the scalar field perturbation 
in the flat slicing from the solutions of the background equations. Very recently the same formula has 
been obtained independently by Kodama and HamazakiB by a different approach. In §§2.3, we use this 
formula to evaluate the curvature perturbation on the comoving hypersurface.-at the end of inflation. 
Then we take the slow rolling limit and recover the result of Sasaki and Stewart a. In §3, we consider an 
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extension of our formalism beyond the limitation of the standard linear perturbation theory. There we 
only assume the smallness of the spatial metric perturbation and the non scalar-type perturbation but 
the other components of perturbations are not supposed to be small. We find almost everything goes 
parallel to the linear case. The summary of this paper is given in §4. 

Throughout this paper, we follow the notation and sign convention used in Kodama and Sasaki u. 

§2. Linear Perturbation 

2.1. Basic equations 

First we write down a set of well known equations for cosmological perturbations. We consider only 
the scalar-type metric perturbation and write the perturbed metric as 

ds 2 = a 2 [-(1 + 2AY)drj 2 - 2B Y i dr^dx J + ((1 + 2H L Y)6 zj + 2H T Y ij )dx i dx j ] , (2-1) 

where Y is the spatial scalar harmonic with the eigenvalue k 2 , Yj = —Je 'VjY, and Yij = k~ 2 V jV jY + 
^SijY. We note the background equations (1-2) and (1-3) can be rewritten in terms of the conformal 
time, n, as 

<^" + 2ft0 p ' + a 2 [/ |p = O, (2-2) 
H 2 = 1 -^' 2 + a 2 U y (2-3) 

where the prime represents a derivative with respect to r/ and TL ■= a' /a. 
The perturbation of the scalar field equation becomes 

X p " + 2H X P ' + k 2 X p + a 2 U lp lqX q - 2{<j) p " + 2H<j) p '}A - cj) p 'A' + <j) p '(3K' - ka g ) = 0, (2-4) 

where x = $4*1 an d we have introduced ka g := H' T — kB and JZ := Hl + \Ht] the former represents the 
shear of the 77 = constant hypersurface and the latter the spatial curvature perturbation. 
The (o)-component of the perturbed Einstein equations is given by 

2 [m 2 A - H(31l' - ka g ) - k 2 TZ] = -</>' • X ' - a 2 U ]p X p + A<j)' 2 - (2-5) 
The (^)-component is given by 

2 [HA - K'} = 0' • x- (2-6) 
The trace part of the Q)-component becomes 

^-(A + TZ)) =4>'-x'-a 2 U\ pX p . (2-7) 







T" + m 


(tt'-i*-.) 


drj 





2 \ a 2 UA + HA' - 
Finally, the traceless part of the Q )-component is 
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ka'g + 2Hka g - k 2 (A + TZ) = 0. (2-8) 

The long wavelength limit of the perturbation, i.e., the limit k 2 /Ti 2 — > 0, is described by taking k 2 — > 
in the above equations. 

An important geometrical quantity which plays a central role in the following discussion is the 
e-folding number of cosmic expansion. The perturbed e-folding number of expansion is defined by 

N= I Hdr, (2-9) 

J T 

where H is the perturbed Hubble parameter given by the 1/3 of the expansion 8 of the n = constant 
hypersurface; 

H = ^6 = H(l + K g Y); 

K g := -A +^(H' L + hB) = -A+ - ha g ), (2-10) 
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and t is the proper time along the integral curve of the vector normal to the r\ = constant hypersurfaces; 

dr = (1 + AY)adn. (2-11) 

One then readily sees that 

/ 1 

(242) 



N = 



H + (TZ? — ^ka g )YJ du. 



Thus the e-folding number will be unperturbed if we take the gauge H' L = B = 0, which implies 
TV = kug/Z. One can then expect that the equations obtained by perturbing the background equations 
by taking N = In a as the time coordinate will closely resemble the perturbation equations in this gauge. 

2.2. N as a time coordinate: the H' L = B = gauge 

If we take N as the time coordinate, the background equations can be written as 



H —r H 4>N 



d 

d~N 
H 2 



3H 2 (j> p N + U^ p = 0, 



1 



:4> 



JY 



~-.u, 



(2-13) 
(244) 



where 4> N = We note that Eq. (244) is the (^-component of the background Einstein equations. 

Together with this equation, the trace part of the (J)-component of the Einstein equations gives 



1 dH 
HdN 



1 2 

77— = -o<Pw> 



(2-15) 



which may be also obtained by substituting Eq. (2T3) into the TV-derivative of Eq. (2T4). 

Let us consider the perturbation equations in the long wavelength limit in the B = H' L = gauge. 
As we noted, in terms of the geometrical variables 1Z and ka g , this gauge condition implies 



W - -ha 



0. 



(246) 



Then Eq. ( |2-4j) reduces to 



H— (H X P N ) + 3H 2 x p N + U^ lqX q + 2U^A - H 2 <j? N A N = 0. 



(2-17) 



Thus among the metric variables, the perturbed held equation contains only A. We also note that A in 
this gauge represents the perturbation of the Hubble parameter, S H/H = —AY, as seen from Eq. (2-10). 

From the ([^-component of the perturbed Einstein equations (2-5), one can see that A is expressed 
in terms of x as 



2UA = —H<p N ■ Xn - U ]P X P , 



(2-18) 



where we have used Eq. (2-3) to simplify the expression. At this point, one may substitute Eq. (248) into 
Eq. (247) to obtain a closed second order equation for X- It is to be noted that only the (^-component 
is necessary to derive the closed equation for x- 



Consider now the trace part of the Q)-component, Eq. (2-7). It reduces to 



2A 



N 



2U{(t>) 
H 2 



A 



U\ pX p 



0~ \p 
n-Xn 



H 2 



(249) 



Together with Eq. ( [2-18 ) , this gives 

A N = d> N -x N - (2-20) 

This equatio n is not independent of Eqs. ( 247 ) and ( |248| ) because of the c ontra cted Bianchi identities. 
In fact, Eq. (2-20) can be directly verified by taking the TV-derivative of Eq. (248) and using Eqs. (2T3), 
(2T4) and (2-17). An equivalent, and perhaps a simpler way to obtain the closed equation for x is to 
eliminate A and A N from Eq. ( |247D by using Eqs. ( |248D and fl2-20p . 

An important fact is that the closed equation for % is obtained without using the (^)-component 
or the traceless part of the (j)-component of the Einstein equations, both of which are absent in the 
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background equations. We also see that Eq. ( |2-20 ) exactly corresponds to the perturbation of Eq. ( 2T5|) . 
This indicates that we are o n the right track. If we obtain a complete set of solutions of Eq. ( 2-17 ) 
(supplemented by Eqs. ( 2-18 ) and (|2-20|) ), the only remaining task is to solve for 1Z and ka g . Let us 
denote the general solution of Eq. ( [2-17| ) by 



(2-21) 



where a runs from 1 to 2n, X( a ) are the In independent solutions and c a are arbitrary constants. 



From the long wavelength limit of the traceless Q)-component, Eq. (2-£), we have 

ka q cx -—. 

a 2 



Under the present gauge condition ( 2-16 ), this implies 

K N = 



1 



(2-22) 



(2-23) 



3aH " a 3 H 

On the other hand, the (°)-component of the perturbed Einstein equations, Eq. ( |2-6| ), is rewritten as 

1 



n 



N 



A - ' X- 



Using Eqs. (2-13), (2-14) and (|2-lq ), TZn can be expressed in terms of x as 

H 2 fdcj)? 



n N = 



2U V dN 



<)>N ■ — 



dx 

dN 



Hence from Eq. (2-23), we find 



W[ X ] ■= 



a 3 H 3 fd<p 



2U 



dN 



■ X 



dN 



N 



(2-24) 



(2-25) 



(2-26) 



should be a constant, which can be directly proved by using Eq. (2-17) with the aid of Eq. (2-18). 

Now, we can calculate 1Z and ka g for each independent solution of the scalar field perturbation. We 
obtain 



N 



N h 



dN 



ka g{a) = 



3W, 



(a) 



(2-27) 
(2-28) 



where Wi a ) is a constant given by 

w (a) = W\x (a) ], 

and Ni, can be arbitrary chosen. We leave it unspecified here. 

The arbitrariness of the choice of Nb is the reflection of a residual gauge degree of freedom in the 
present gauge. Note that there is a trivial solution x — for which 7Z = c° = constant. The addition of 
this solution corresponds to a variation of Nb. Hence we have 2n + 1 integration constants for the whole 
set of the perturbation equations in this gauge. This is not a contradiction because there is a residual 



gauge degree of freedom in the gauge B 
coordinate, 

Hl transforms as 



H' L = 0. Under an infinitesimal transformation of the time 
N N - SN, (2-29) 



Hl —* Hl + SN. 



(2-30) 



Hence the gauge condition H' L = allows a further gauge transformation given by SN — c = constant, 
which corresponds to an infinitesimal time translation mode. Applying this gauge transformation to the 
null perturbation, we obtain a pure gauge mode, 



X = cc/) N , TZ = c. 



(2-31) 



This fact implies that one of the 2n solutions of Eq. (2-17) should be proportional to this time translation 
mode. In fact, we can verify that % = <p N is a solution of Eq. (2T7) by direct substitution. Hence we 
may set 

= <V (2-32) 
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Then we also find Wm = 0. 

The issue of the number of physical degrees of freedom becomes transparent by constructing a gauge 
invariant quantity. A convenient choice is the scalar field perturbation Xf on the flat hypersurface 
(defined by 1Z = 0). It is given by 

Xf = X - <P N K. (2-33) 
Then the general solution for Xf 1S given by 

Xf = ? a (x (a) ~ W {a) il>) , (2-34) 

where = and we have redefined c 1 by the replacement; c 1 — c° — ► c , and 



N 



dN 



t^^kWH' (2 ' 35) 

One sees that the change of Nb is always absorbed in the redefinition of c 1 . Hence, only the 2n integration 
constants remain. 

Since there exists at least one solution that has a non- vanishing W, we may assume ^ without 
any loss of generality. Then it is worthwhile to mention that it is always possible to set = for 

a ^ 2 by redefinition of the complete set of solutions. In fact, W vanishes for the linear combination of 

the solutions defined by X( a ) ~ %(2) f° r a — 3- ^ mav ^ e a ^ so useful to note that c a W( a ) describes 

the amplitude of the adiabatic decaying models, while the adiabatic growing mode amplitude is given 
bye 1 . 

Let us now consider the perturbation of the background equations (2-13) and (2T4). Let us assume 
that the general solution for the background equations is known. Except for the trivial time translation, 
the general solution contains 2n — 1 integration constants, A a , which distinguish the various solutions. 
Each solution gives a curve in the phase space of the scalar field parametrized by some time coordinate 
A 1 . Let x A — (x P i 7r9 ) [A — lj 2, * - • ,2n) be the phase space coordinates where ir q are the momentum 
variables. Then the general solution is expressed as x A — x j4 (A' 5 ). We may regard this as a coordinate 
transformation of the phase space coordinates. We see that the perturbation of a given background 
solution is given by the Jacobian dx A /d\@, including the time translation mode. It is then easy to 
convince ourselves that a necessary condition for these background solutions to describe the solutions 
of the perturbation in the long wavelength limit is the non-disturbance of the time coordinate under 
the presence of a perturbation, since the commutable property of the partial derivatives; d 2 / 'd\ a d\@ = 
d 2 /dX^dX 01 , should be maintained. As we have seen, the e-folding number N = In a is indeed such a 
coordinate for the gauge B = H' L = 0. Thus we take iV as the time coordinate in the background 
equations and set A Q = (N, X a ). Note that the background equations (2T3) and (2T4) have no explicit 
iV-dependence. Hence <p N is a solution to the perturbed background equations, which corresponds to the 
time translation mode. 

Taking the derivative of Eqs. (2T3) and (2T4) with respect to A", we obtain 



2U JT = H2<t>N ' Jn^ + {2 ' 37) 

where the suffix A represents d/dX a . We find these equations are equivalent to those for the perturbation 



in the long wavelength limit in the B = H' L = gauge, Eqs. (2T7) and (2TE), respectively, with the 



identifications 4>\ = X an d H\/H = —A. Further, as we have already seen, the A a -derivative of Eq. (2-15) 



is equivalent to Eq. (2-20). Thus we conclude that a complete set of solutions in the long wavelength 
limit, \{a)i can be constructed from the solutions of the background equations; X( a ) = 94>/dX a . Once 
X( a ) are obtained, the corresponding set of gauge-invariant quantities Xf(o) IS readily obtained as 

Xf{c) = X( a ) ~ W ( a)i>, (2-38) 

where W( a ) is given by Eq. ( 2-28| ). 
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Before closing this subsection, we mention the g enera lization of the scalar field metric h pq to a non- 



trivial one. Except for the perturbed field equation (|2T7| ) or (2-36), the only modification is to replace 
all the derivatives in the equations with the covariant ones with respect to h pq . For example, 

where r p qr is the connection of h pq . As for the perturbed field equation, we have to add the curvature 



term in addition to the covariantization of the derivatives. That is, Eq. (2T7) is modified as 



(H X P N ) + 3H 2 x p N + H 2 RV rqs ^ s N x q + U^ lqX q + 2U^A - H 2 <j? N A N = 0, (2-40) 
and similarly for Eq. (2-36). 

2.3. Curvature perturbation on the comoving hypersurface 

Among the various geometrical quantities, one of the most convenient representations of the per- 
turbation amplitude is the curvature perturbation of the comoving hypersurface, _K. C . For the adiabatic 
growing mode perturbation, it is known to stay constant on super-horizon scalesEP. Here, we relate 1Z C 
with the scalar field perturbation and present a formula that can be used to evaluate 1Z C at the end 
of inflation in terms of the initial data of the scalar field perturbation. Then we take the slow rolling 
limit and show how our formula reduces to the one obtained in Ref. |]). To obtain 1Z C that is directly 
relevant to observational quantities, one must solve the evolution of 1Z C during the reheating stage after 
inflation. One may also have to evaluate the possible contribution of the isocurvature perturbation. We 
defer these issues to future studies and derive here a formula valid up to a time N = N e before the 
reheating commences. It may be noted, however, that there are situations in which one can neglect the 
time variation of 1Z C during reheating. For example, if all the trajectories in the phase space of the scalar 
field converge to one path before reheating, the evolution during reheating is essentially the same as the 
case of single-scalar inflation and there will be no additional change in the amplitude of 1Z C . 

As we have seen in the previous section, the spacetime configuration of the scalar field on super- 
horizon scales is completely determined by the background solutions and the parameters A Q characterizing 
the background solutions play the role of phase space coordinates. This implies the value of the scalar field 
depends on the spacetime coordinates only through the phase space coordinates for the spatially 

homogeneous configurations; <p{x^) — <p(X a (x^)). Thus it will be useful to express 1Z C in the language 
of the phase space of the background scalar field. 

The comoving hypersurface is determined by the condition, = 0, i.e., 

0jv • dcf> = 0, (2-41) 

In general, we cannot expect this condition to determine a surface in the phase space of the scalar field 
since it is not integrable. However, in the case of linear perturbation, this condition reduces to 

4> N (N)-(<f>-j>(N)) = 0, (2-42) 

where the barred quantities represent the background values. Now it is manifest that this condition 
determines a surface in the phase space. For each A a , the value of N at which the comoving surface in 
the phase space is crossed is different. We denote this difference by AN. Namely, 

4> N (N) ■ (0(JV + AN, X a ) ~ 4>(N)) = 0, (2-43) 

where AN is a function of A" = (N, X a ). From the spacetime point of view, through the dependence of 
A Q on x^ , AN is a function of x^ . Hence we may consider the infinitesimal coordinate transformation 
N — > N — AN and move from the H' L = B = gauge to the comoving gauge. Then we have 

K C (N) =n(N) +AN. (2-44) 

On the other hand, the comoving condition ( 2-43| ) is rewritten as 



4> N (N) ■ ( X (N) + 4> N (N)AN) = 0. (2-45) 
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Thus we have AN 



'{4>n ' X)/4>n- Therefore we obtain 



c a W, 



(a) 



N 



dN 



:SiV 



(2-46) 



where we remind that X( a ) = d(j)/d\ a . Given the background solutions, this gives the evolution of 1Z C 
up to N = N e . 

In the above, we have derived the expression for 1Z C in the language of the phase space for the 
background solutions. But, of course, it can be derived in the language of the standard linear perturbation 
theory. In the defining equation for Xf^ Eq. (2-33), we evaluate the right hand side on the comoving 
hypersurface to obtain 

Xf = Xc-^n^c (2-47) 

where x c i s the scalar field perturbation on the comoving hypersurface. Then taking the inner product 
of the above equation with <f> N and using the comoving hypersurface condition <f> N • X c = 0, we find 



n r = - 



• Xf 

4> 2 n 



(2-48) 



Inserting the expression ( |2-34 ) for Xf to this equation gives Eq. (p46|) . 

The remaining task is to determine c a . In the inflationary universe scenario, scalar field perturbations 
are generally induced by quantum vacuum fluctuations. The evaluation of the vacuum fluctuations is 
most conveniently done with respect to Xf H n ce there exists a closed action for Xf that gives a complete 
description of the scalar- type perturbation EF. Hence we assume that Xf ^ s given at some initial epoch 
when the wavelength exceeds the horizon scale. 

Suppose that the initial condition of Xf an d 
can be readily evaluated by inserting Xf m t° the operator W: 



DXf 

— is given at a time, N = Nn. The constant c a W/ n ) 

dN w 



Hence, 



c a W, 



(a) 



a 3 H (D<j> 



N 



6 \ dN 



Xf - <f> 



N 



3iP 

DXf 

dN 

A-t. 



c a W (a) 



Then ado pting the notation x := \ for ^4 < n and \ := \n ™ for ^4 > 
Eq. (HI) as 



(X- 1 )^ {xf + (c f3 W (p) )^} 



N 



(2-49) 

(2-50) 
is evaluated from 
(2-51) 



where Xf an d ip ar< 3 defined in the same way as x, an d (x 1 )a ^ s the inverse matrix of xf a )i ^ e 



(2-52) 



So far our discussion has been completely general. Now we consider the slow rolling limit. We assume 
all the components of the scalar field satisfy the slow rolling condition: 



4>N « 



In this limit, the background equations reduce to 

jj\p 



D4> 



N 



3iF 



dN 



(In 17) I 



« 10 



N 



H 2 



-U. 



(2-53) 



(2-54) 



Thus the momentum variables cease to be dynamical and the 2n-dimensional phase space reduces to the 
n-dimensional configuration space of the scalar field. Correspondingly the number of the parameters A" 
reduces to n. The perturbed field equation is very simple: 



dN Vx 



05 = -On £7)1 V' 



A ' 



(2-55) 
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where we have assumed the scalar space curvature is small; |i? p r <;s| <C 1.0 From this equation, we 
easily see that VF[</>>J = for all the n-independent solutions. Thus the slow rolling assumption kills the 
adiabatic decaying mode as in the case of single-scalar inflation. This implies that \ F is directly given 
by the background solutions: 



Xf = c X( a ) 



dX a 



(2-56) 



where a now runs from 1 to n. In other words, the flat slicing does not perturb the e-folding number N 
in the slow rolling limit. 

The formula ( |2-51 ) for c a reduces to 



N 



dX a 



P 

Xf 



(2-57) 



N„ 



For definiteness, we set X(i) = 4>n (i- e -> A 1 = N) as before. The curvature perturbation 1Z C is then given 

by 



a 4>N ' X(a) 



N 





'ON / 




'd\ a p 


&N ■ 4>,a \ 




d(/)P Xf _ 


+ 


d4>p Xf _ 


N a 4>N J 




N 



(2-58) 



where a = 2, 3, ■ • • ,n and <fi a = d<p/d\ a . 

If we synchronize the time to be N — N e at the end of inflation for all the background scalar field 
trajectories, we have dU/d\ a = at N = N e . Then using the slow rolling equation of motion (2-54) we 
find 

(2 ' 59) 

at N = N e . Thus the evaluation of Eq. (2-58) at N = N e gives 

~dN 



n c {N e 



d<j)P 



P 

X F 



(2-60) 



No 



This formula coincides with the one derived by Sasaki and Stewart t 



§3. Quasi-Nonlinear Perturbation 

An advantage of the formalism developed in the previous section is that it can be extended so as 
to take into account the non-linearity of the scalar field dynamics. The basic idea is to consider the 
linearization with respect only to inhomogeneities of the spatial metric. For simplicity, we again assume 
the trivial metric for the scalar field space, h pq — A 
We consider the following metric form, 



Jpq- 



ds 2 = --L;dN 2 + e 2N ({l + 2H L )Sij + 2H Tij )dx l dx\ 



(3-1) 



where H depends both on N and x 1 . Here and in what follows the tilde is attached to a non-linear ly 
perturbed quantity. We consider and -ffrij to be small of O(e). In accordance with the linear case, 
we assume the existence of a potential function Ht for Hxij such that 



3 J 



Hi 



(3-2) 



*' It seems this condition is necessary for the slow rolling approximation to be consistent. However, we are 
unable to give a proof. 

**) In Ref. the minus sign on the right hand side of this formula is absent. The reason is that N there 
is defined as the e-folding number counted backwards from the time N — N e for convenience of the spectrum 
evaluation. 
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We also assume the spatial derivatives of H and <p are small of 0(e). This assumption corresponds to 
taking the long wavelength limit. Furthermore, we assume the gauge condition dH^/dN = 0(e 2 ) can be 
imposed consistently. We then linearize the scalar field and Einstein equations with respect to e while 
keeping the other nonlinear terms as they are. 

First let us consider the scalar field equation. Keeping the terms linear in e, the field equation reduces 

to 







1 



--d„ 



9" 



-gd„4? 



JL (1 _ 3A) ^» (1 + 3ffl) _i> 



4> p + u{4>yv 



(3-3) 



This is the same form as the background equation (2T3). 

Next let us consider the (g)-component of the Einstein equations. To the linear order in e, we have 



Lr o 



-3H 2 , 

T\ = ~H 2 ~^ N -UQ>). 



Therefore we find 



H 2 1 



~ & 4>N 



(3-4) 
(3-5) 



which is also the same form as the background equation (2T4). 

Thus we find that, in the long wavelength limit, solving the equations for the perturbed scalar field 
in the dH^/dN = gauge is completely equivalent to finding the spatially homogeneous background 
solutions. 

The remaining task is to determine 1Z := Hi + \Ht and ka g :— e N H dH^ / dN as before (note that 
k<j g is not k times <j g but regarded as a single symbol here). The traceless part of the (J)-component of 
the Einstein equations is essentially the same as Eq. (2-8) also in the present case. This tells us that we 
have 



1 



ka g = 3e HTZn oc , x . 
The (^)-component of the Einstein tensor becomes 

G°i = ViH 2 + 2H 2 ViK N 
while that of the energy momentum tensor is given by 

T°i = -H 2 cj) N ■ Vi0. 

Hence we have 

^HV.TZn = -e 3N (v>H +^H4> N - Y ; cb 



(3-6) 

(37) 

(3-8) 
(3-9) 



which is expected to be constant in time because of Eq. (3-6). In fact, we find 
d 



„3iV 



dN 



Hv l n N ) 



d 
dN 

3N I 



„3N 



e 
2H 
0. 



2(7(0) 



-Hcf) N ■ V 4 



H 



ViH + H'ct> N V l( f> N + ViU(4>) 



(3-10) 



where Eqs. (3-3) and (|3^) are used in the second equality and the spatial derivative of Eq. (|3-5| ) in 
the last equality. This verifies the consistency of our assumptions, in particular the gauge condition 
dH L /dN = 0(e 2 ). 

Let us reinterpret the discussions given in §§2.3 in the present context. For simplicity we assume 
that the linear perturbation is valid at N — Nq. Otherwise, the evaluation of the initial perturbation 
becomes too difficult. As in §§2.3, we give the initial data in terms of Xf an d dxp/dN. 
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Then from Eq. (3-6), the evolution of 1Z is given by 



n = c n 



N 



dN 



I No e^H 1 

where the coefficient C-ji is expressed in terms of the initial values of Xf anc ^ dxp/dN as 



C-r. — 



e iN H (d<t> 



(i 



dN 



■ Xf - <P 



N ' 



dx F 
dN 



(3-11) 



(3-12) 



J JV 



which corresponds to Eq. ( |2-50 ) for the linear case. In the above, we have chosen TZ(Nq) = by setting 
A*6 = iVo for co nveni ence. 

Using Eq. ( |2-34 ) , the initial data for the scalar field are given by 



4>{N ) = 



<>(No) + Xf(No), 
Ck 



1 



e 3N °H 



&n(N )+Xfn(No), 



(3-13) 



where Ho = H(Nq) and X^^v = Axf/^N. The evoluti on of <fi is determined by solving the nonlinear 
background field equation (3-3) supplemented with Eq. ( |3-5| ). 

In §§2.3, to evaluate the amplitude of the perturbation at the end of inflation, it was necessary to 
move to the comoving hypersurface. In terms of the background solutions, this amounts to finding a 
comoving surface in the phase space. However, as noted there, the comoving condition (2-41) does not 
specify a surface in the phase space in the nonlinear case. To circumvent this difficulty, here we propo se to 
use the constant Hubble hypersurface in substitution of the comoving one. It is apparent from Eq. (3-5) 
that there exists a surface in the phase space corresponding to the constant Hubble hypersurface. From 
the fact that the right hand side of Eq. (|3-9[) is constant in time, we find 



(3-14) 



VilnH = -~d> N ■ Vi<j>-\ 

1 2 e 3N H 



where dj is a time-independent vector. At a later epoch, the second term in the right hand side can 
be neglected and the difference between the comoving gauge and the constant Hubble gauge becomes 
negligibly small. 

Then as a function of Xf(^o) and Xfn(No), we define ANh as the difference of the time to cross 
the constant Hubble surface in the phase space, namely, by the condition, 

H (cj)(N + AN H ), (t> N (N + ANh)) = independent of A a . (3-15) 

Finally, we find the curvature perturbation on the constant Hubble hypersurface TZh as 

n c ~n H = n + AN H - f>i6) 

In the case the slow rolling condition ( 2-5 3| ) is satisfied at the initial time, the first term in the right 
hand side of Eq. (3T6) can be totally neglected and we simply have 



Tl c ~ TZ H = AN H , 
where ANh will now be a function of Xf(Nq) alone. 



(3-17) 



§4. Summary 



In this paper, we have investigated the dynamics of a multi-component scalar field on super-horizon 
scales, i.e., in the long wavelength limit. We have shown that there is a simple relation between the 
perturbation equations in the long wavelength limit and the background equations. That is, the derivative 
of the general solution of the background equations with respect to a parameter that characterizes different 
solutions satisfies the same equation as the perturbation in the long wavelength limit does. However, we 
have also found that the explicit form of the relation depends on the choice of gauge, and the choice of 
gauge corresponds to that of a time coordinate in the phase space of the background scalar field. We have 
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found that the simplest form of the relation is obtained in the gauge in which the e-folding number N of 
cosmic expansion is unperturbed even under the presence of the perturbation. Then using this result, we 
have given a method to calculate the amplitude of the spatial curvature perturbation on the comoving 
hypersurface, 1Z C , from the knowledge of the background solutions alone. 

As a natural extension of our approach, we have considered to take into account the nonlinearity of 
the scalar field dynamics in the perturbation. We have found that this can be actually done under several 
reasonable assumptions. The result provides a powerful tool to evaluate the effect of nonlinearity of the 
scalar field potential during the inflationary stage. In particular, the effect of non-gaussian statistics 
of the perturbation, given a gaussian distribution of the initial perturbation due to quantum vacuum 
fluctuations, can be evaluated by studying the background solutions alone. n 

The present nonlinear -extension is similar to the so-called anti-Newtonian approximation u or the 
gradient expansion methodlMP. It may be worthwhile to clarify how our result is related to these methods. 
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Appendix A 

Other Choices of Gauge 

In this appendix, we compare the long wavelength perturbation equations and the background equa- 
tions in a couple of other choices of gauge. We find the choice of the time coordinate in the scalar field 
phase space determines the corresponding gauge condition for the perturbation as in the case of the 
H' L = B = gauge. However, the resulting equations turn out to be less tractable. 

A.l. t as a time coordinate: the synchronous gauge 

In the synchronous gauge, A = B = 0, the cosmological time t remains to be the proper time along 
curves normal to the constant time hypersurfaces. Hence we expect t to be the relevant time coordinate 
to parametrize the background solutions. 



The perturbation equations (2-4) and (2-5) in the long wavelength limit in the synchronous gauge 
become 

X P + 3Hx p -^{^- + ulP WX q = 0. (A-1) 
2H Un - £ZL\ = • ± + U lpX p ■ (A-2) 



/ • fc<7„ \ 

Eliminating 37?. from these equations, we obtain an equation written in terms of x alone. 

V a J 

Note that ( 1Z — ) in this gauge represents the perturbation of the Hubble parameter as seen from 

V 3a J 

Eq. (2-10). 

Once we know the solutions for we can evaluate the spatial curvature perturbation with the aid 
of the (^)-component of the Einstein equations, Eq. (2-6), which now reduces to 

K = -^-X. (A-3) 



Integrating this equation, 1Z can be evaluated. But different from the H' L = B = gauge, we have no 
knowledge of the explicit time dependence of the integrand at all. So the integration must be done for 
each mode separately. 

On the other hand, the time dependence of ka g can be found from the traceless part of the (*■)- 
component, Eq. ((2~§1), as 

kdg cx . (A-4) 
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Using Eqs. (A-2) and (A3), this implies 

a 3 ■■ 

kv g a 2 = 2Tj(<f> • X ~ 4> • X) = constant . (A-5) 

As before, the general solution of Eq. (A-l) with (A-2) contains 2n integration constants. When we 
evaluate 1Z, there appears an additional integration constant. Hence, we have 2n+l integration constants 
in total. This situation is exactly parallel to the case of the H' L = B = gauge. Namely, one mode is 
responsible for a gauge degree of freedom, since the synchronous gauge condition allows an additional 
gauge transformation given by t — > t + c, where c is a constant. Then the gauge mode is found to be 

X = c0, 11 = cH. (A-6) 

Also parallel to the H' L = B = gauge, there exists a trivial solution: 

X = 0, TZ = constant. (A- 7) 

As before, this gives one of the solutions for the gauge-invariant scalar field perturbation \f as 

Xf = X~ ^n<x<t> N . (A-8) 

Now let us turn to the comparison of the perturbation equations with the background equations. 
As we have mentioned in the beginning, since the cosmological time is unperturbed in the synchronous 
gauge, the relevant time coordinate will be t. Thus taking the derivatives of Eqs. (1-2) and (1-3) with 
respect to A" = (t, A a ), we have 

4>{ + 3H x cj>P + 3H<P P X + U\ p \ q <t> q x = , (A-9) 

and 

2HH X = i (<t> ■ A + U [P 4>\) , (A-10) 

where the suffix A again represents the partial derivative with respect to A Q . It is readily seen that these 
are equivalent to Eqs. (A-l) and (A-2), with the identifications, 

A = X, H x = U-\^. (A-ll) 
3 a 

A. 2. H as a time coordinate: the constant Hubble gauge 

As another example, we consider to take the Hubble parameter H as a time coordinate. Then the 
relevant choice of gauge will be to take the constant Hubble gauge in which H is unperturbed under the 
presence of the perturbation. 

From the expression for the perturbed Hubble parameter, Eq. (2-10), we find the condition for the 
constant Hubble gauge as 

HA-W + \k<j g = Q. (A-12) 

In this gauge, the perturbation equations ( |2-4j ) and ( |2-5[ ), in the long wavelength limit, reduce to 

X p " + m X p ' - {W + n^ 1 } A - <f'A' + a 2 U^\ qX q = 0, (A-13) 

A = \ {4>'-x' + a 2 U lp x p ). (A-14) 




From Eq. (2-8), we find ka g cx 1/a 2 , and from Eq. (2-6), we have — (2/3)ka g = cj> Hence, we find 



3a 2 

ka g a 2 = —(f) 1 ■ x = constant. (AT5) 



Taking the time derivative of this equation, we find the constraint, 

(/>' ■ X ' = a% P X p ■ (A-16) 
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This means that the initial condition for the scalar field perturbation cannot be chosen arbitrarily. This 
is because the constant Hubble gauge condition completely fixes the time slicing and hence the constant 
time surfaces in the phase space of the scalar field. So the perturbation in the direction normal to the 
constant time surfaces is not allowed. 

Hence, in the present gauge, the general solution for scalar field perturbations has 2n — 1 integration 
constants. The curvature perturbation can be evaluated by integrating 



n' = HA-ld>'-x = 2n ^f -\<t>''X- 



(A-17) 



Again, the integration of this equation is non-trivial. 

As in the previous two cases, we have a trivial solution 



X = 0, 1Z = constant. 



(A-18) 



This again tells us that cj) N is a solution for Xf- 

Now let us consider the background equations. Using H as the time coordinate, the background 
equations are written as 



1 d <fi 



t 3l_ & JL^ 

</>H dH <f>H 0ff " 



H 2 



4> 2 H 



u 



U^ p = 0, 



(A-19) 



where subscript H represents the derivative with respect to H. The latter equation is the equation that 
constrains the scalar field in its phase space. Different from the previous two cases, these equations do 
not have the invariance with respect to the time translation because they contain H explicitly. Hence, 
(p H is not a solution of the perturbation equations. 

Taking the derivative of Eqs. (AT9) with respect to A = A a , we obtain 



2H4^ - {2<t?" + H<jf'} 



0' ■ 0A = a 2 U\ p (, 



20' ■ <t>'\ 
0' 2 



20' ■ 0'a 
. 0' 2 



a 2 C/V A = 0, 



(A-20) 



These are equivalent to Eqs. (A- 13) and (A- 14) supplemented by the constraint (A- 16) 
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